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The question addressed is whether magnetic materials based on physical small world networks 
are possible. Physical constraints, such as uniform bond length and embedding in three dimensions, 
are the new features added to make small world networks physical. Results are presented to further 
determine if physical small world networks can exist, and the effect of the small world connections 
on the critical phenomena of Ising models on such networks. Spectra of the Laplacian on randomly- 
collapsed bead-chain networks are studied. The scaling function for the order parameter of an Ising 
model with physical small world connections is presented. 



I. INTRODUCTION 

It is well known that novel magnetic behavior is 
present in systems governed by fixed points in other 
than three dimensions, such as quasi-one-dimensional 
and quasi-two-dimensional magnetic systemsi Recently 
there has been renewed interest in the behavior of com- 
plex networks, 2 including small-world (SW) networks. 3 
These networks mainly originate in areas in technology, 
or the social sciences^ Properties of networks have been 
used in applications, for example, to devise perfectly 
scalable parallel algorithms. 5 Magnetic models, includ- 
ing i s i n gSigi2ii2iiiiiS models, have been studied on such 
networks. Recently, the question has been asked whether 
or not materials based on physical SW networks are 
possible^ Physical constraints, such as uniform bond 
lengths and node sizes, and that the network must be 
embedded in three dimensions, are the new features we 
add to SW networks to make them physical?! 3 . Since mag- 
netic models on such SW networks are usually governed 
by a (slightly modified) mean- field fixed point ji 2 ? simi- 
lar to mean-field fixed points for polymer systemSfi^i 5 ^^ 
magnetic spins on physical SW networks should exhibit 
novel magnetic behavior. 

Here we present results to further determine whether 
physical SW networks can exist, and the effect of the 
SW connections on the critical phenomena of Ising mod- 
els on such networks. We present results of a project 
with randomly-collapsed bead-chains to determine the 
topological properties and the properties of the Lapla- 
cian on such networks. The Laplacian operator is the 
most fundamental operator governing physical behavior 
on the network, such as diffusion or collective excita- 
tions in various interacting systemspiii^ The starting lat- 
tice for the Ising simulations are either one-dimensional 
chains or two-dimensional toruses, with the physical SW 
connections added to the underlying lattice. We show, 
for example, how the mean-field behavior of Ising fer- 
romagnets arises as the lattice size and density of SW 
connections is varied. 



II. MAGNETIC MODELS 

The magnetic models simulated have Ising spins, 
5j=±l, placed on the sites of a one (d=l) or two di- 
mensional (d=2) square lattice, with periodic boundary 
conditions and nearest neighbor ferromagnetic exchange 
coupling J\. A selected number of SW bonds, with fer- 
romagnetic coupling J2=4Ji, are also included (Fig. 1). 
The Hamiltonian is 

"H = —J\ ^ ' SjSj — J% Saw jS sv/ j. (1) 

Our Monte Carlo simulations^ use a Glauber dynamic 
with random, single-site updates. When there are no SW 
bonds for the d=2 model the critical temperature ksT c « 
2.269Ji. A modified Wang-Landau method 2 — has been 
used to obtain comparisons with standard Monte Carlo 
results. 

III. PHYSICAL MODELS 

Two methods to construct physical networks, close to 
SW networks, were used. In both cases, the networks 
formed a physical network because the network is em- 
bedded in d=3 in such a way that all bond lengths are 
short-ranged and nearly equal. 

In the first method, four high school students sequen- 
tially numbered all beads, randomly collapsed, and glued 
Mardi Gras necklaces of two different lengths. The bead 
size was 8 mm, and the lengths of the necklaces were 
106 cm (about 60 beads) and 254 cm (about 155 beads). 
Beads were cut from the necklaces, with a record kept 
of the adjacency matrix, A. The matrix A has 1 as its 
i,j element if bead i is connected (either by the neck- 
lace string or by glue) to bead j. Otherwise the matrix 
element is zero. 

Physical SW networks were also constructed by imple- 
menting a simulated annealing program to find an en- 
ergy minimum of an arrangement of atoms starting with 



FIG. 1: Examples of the lattices with SW connections that were studied. The filled circles are the locations of the Ising spins, 
the light solid bonds have strength Ji and the dashed bonds (the SW bonds) have strength J2. (a) A one-dimensional SW 
graph with N = 20 and five randomly chosen SW bonds, (b) A square-lattice with L=4 with L SW bonds. 



a lattice such as in Fig. 1(a). The energies included were 
those due to the bond lengths, bond angles, dihedral an- 
gles, a Lennard- Jones interaction, and partial ionization. 
Energy parameters were taken from the AMBER? 1 data 
set. These physical SW networks will be analyzed else- 
where, but lead to the question of how the number of SW 
bonds scale with the system size. 

IV. DATA AND ANALYSIS 

The ordered (non-zero) eigenvalues of the (negative) 
Laplacian on the glued Mardi Gras bead-chain networks, 
— Aij + 5ij^2 £ Ai£ [Fig. 2], were found to be very sim- 
ilar to those of SW networks, such as in Fig. 1(a). In 
particular, the value of the smallest eigenvalues (govern- 
ing, e.g., low-energy excitations of such systems) are sig- 
nificantly larger than those in the spectrum of a pure 
Id necklace (where with N sites the eigenvalues are 
2 [1 - cos (27rfc/A0] with k = 0, 1, 2, • • • , N - 1). This 
indicates the opening of a pseudo-gap 17 ' 18 and the corre- 
sponding mean-field-like behavior in the physics of fused 
Id systems, commonly observed in SW networks^ 

For physical SW networks, one unanswered question 
is how the number of SW bonds scales with the system 
size. Consequently, we have investigated cases in which 
a finite lattice has a reasonable number of SW bonds, 
but the density of SW bonds goes to zero as the system 
size increases. The order parameter m should scaled as 
to = L-P' v F{y) with y = L x l v t = L 1 /" \(T - T c )/T c \. 
For large y with L — > 00 and t <gc 1 the scaling function 
should asymptotically be F(y) w By 13 for T < T c and 
F(y) w B'yP~ v for T > T c . Since the density of SW 
bonds goes to zero as L — > 00, the critical exponents and 
critical temperature should be those of the underlying 



Ising model. Fig. 3(a) shows the scaling function for to 
for &LxL square lattice with no SW bonds, and Fig. 3(b) 
for L SW bonds of strength J2 = 4Ji. Although asymp- 
totically the agreement with the normal Ising exponents 
is apparent, the SW bonds modify the scaling function. 
This illustrates that some properties of the mean-field 
behavior due to SW bond cfi 12 ' 13 should be observed in 
these finite (hundreds of nanometer in length) systems. 



V. CONCLUSION 

We have presented initial investigations of physical SW 
networks, and the magnetic behavior caused by SW in- 
teractions in Ising systems. Further work is required to 
definitively illustrate the behavior of magnetic systems 
on physical SW networks. Even if physical SW connec- 
tions are shown to be irrelevant in the thermodynamic 
limit, the scaling functions for nanomaterials change sub- 
stantially, providing quasi- SW behavior. Consequently, 
the current known results should provide a strong impe- 
tus to attempt to synthesize (or construct) and measure 
magnetic properties of such systems. 
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FIG. 3: The scaled order parameter for Ising models on L 2 square lattices. The asymptotic slopes have values of ft=\ and 
ft — v=— g. (a) for no SW connections, (b) for L SW connections, as in Fig. 1(b), with strength J2=4Ji. 
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